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Abstract — This paper considers the use of punctured con- 
volutional codes to obtain pragmatic space-time trellis codes 
over block-fading channel. We show that good performance can 
be achieved even when puncturation is adopted and that we 
can still employ the same Viterbi decoder of the convolutional 
mother code by using approximated metrics without increasing 
the complexity of the decoding operations. 

I. Introduction 

A relevant result obtained in wireless communications is 
that the use of Space-Time Codes (STC) with both mul- 
tiple transmitting and receiving antennas can be exploited 
to mitigate the effect of fading without sacrificing spectral 
efficiency. Several research activities in the last decade has 
been devoted to find STC able to improve the performance 
in terms of achieved diversity and coding gain [1]. Since the 
introduction of STC, the problem of determining the best STC 
has been always a difficult task, especially for a large number 
of transmitting antennas. 

In [2] a new approach to space-time coding called Pragmatic 
space-time coding (PSTC) has been devised. The pragmatic 
approach consists in using common convolutional codes as 
STC over a MIMO channel; it has also been shown that PSTC 
could be used with both QPSK and BPSK to achieve maxi- 
mum diversity and good performance. The use of convolu- 
tional codes as STC allows decoding by using the same Viterbi 
decoder with only a proper change in metrics evaluation. 

However, in the case of high-rate codes, R — k/n, the 
number of operations that the decoder has to perform grows 
exponentially with fc, and the number of paths to be stored 
increases rapidly. Puncturation was initially introduced for 
convolutional codes to avoid the complexity issue of the 
Viterbi decoder in case of high rate codes. Puncturation 
consists in deleting one or more bits of a codeword; in this way 
rate pj (p+m) punctured code can be obtained by periodically 
puncturing a low rate 1/n convolutional code, i.e. by erasing 
to bits for each period of length p. 

In this paper, following the same aforementioned pragmatic 
approach to STC, we will introduced puncturation, using 
therefore punctured convolutional codes as ST codes. We will 
call this family of codes punctured pragmatic space-time codes 
(P 2 -STC). The use of puncturation in conjunction with STC 
has also been investigated in [3], where the pragmatic approach 
was not adopted and only the case of two transmitting antennas 
was addressed by designing the codes to preserve full diversity 
gain only for short error sequences. 



II. Pragmatic space-time codes: 

In this work we consider P 2 -STC in a block fading channel 
(BFC)(see [4] [5]) where fading coefficients for each couple 
of transmitting and receiving antennas are constants in blocks 
of B bits and independent block by block. L different blocks 
are experimented by a codeword. In this situation, being N 
the number of transmitting antennas and M the number of 
receiving antennas, we can generalize the Singleton bound for 
BFC to obtain the following diversity bound 



div < 1 + [LN(1 - R)\ 



(1) 



where div is the maximum achievable diversity degree per 
receiving antenna. 

This means that, by puncturing a rate 1/2 convolutional 
code in quasi-static fading channel (L = 1) with N = 2,M = 
1, we cannot achieve div = 2 since R > 1/2. For a rate 1/3 
convolutional code with N = 3, M = 1 and L = 1, we cannot 
achieve div = 3 if R > 1/3. Our results prove however, that 
we can still implement puncturation and preserve div = 2 if 
l/3<i?<2/3. 

To describe a system with P 2 -STC we refer to the block 

scheme in Fig[T| We denote by = [c^ ■■■ c ^n\ T a su per- 
symbol, i.e. a vector of symbols simultaneously transmitted by 
the N transmitting antennas at the time t. A frame is composed 
by a sequence of F super-symbols (every T s seconds, N 
symbols are sent in parallel on the N transmitting antennas). 
If we indicate the complex envelope of the signal received 
by antenna s at time t with r^p and with a® the fading 
coefficient between antenna s and i, at time t, we have: 
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Fig. 1. Block scheme for pragmatic space-time codes (transmitter side), 
4— states, lbps/Hz, BPSK modulation, n = 2 antennas, (5, 7)g generators, 
Decoder stands for Viterbi decoder with proper branch-metrics. 



To perform maximum likelihood decoding, the Viterbi algo- 
rithm should estimate the input bit sequence b in this way: 



A? 



b = arg min £ £ | r W - £ a«cf (6) ^ 



(3) 



where q*' (b) are the transmitted symbols corresponding to the 
path in the trellis labeled by the bit sequence b, and 



JV 



(4) 



are the metric increments. The fundamental problem is, how to 
compute these metric increments when we adopt puncturation. 
Consider for example a rate 1/2 code used with the punctura- 
tion matrix 

with two transmitting antennas and BPSK modulation. 
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Fig. 2. Relationship between received symbols and coded bits for a P 2 -STC 
obtained from a R=l/2 convolutional code with N=2, M=l 



Fig. [2] shows the relationship between the received symbols 
A*' and the coded bits. We can see that the second and 
third transition of the trellis diagram are labeled with only 
one coded bit, i.e., symbol A 2 * 1 carries coded bits from two 
different transitions and therefore we must use the same 
received symbol A 2 ^ to compute two transition metrics on the 

(2) 

trellis. Since the coded bit c\ cannot be recovered from the 
trellis diagram at instant t = 2 we can use it only for the third 
transition metric. To implement ML decoding we should join 
two transitions into one such that no symbol contains bits from 
two transitions anymore, as already pointed out in [6]. This 
however, would mean that we could not use a conventional 
Viterbi decoder of the rate 1/2 code. To exploit the same 
Viterbi decoder we must actually choose an approximation 
of ML decoding that uses approximated metrics. In [6] a 
very similar problem arise in the context of trellis coded 
modulation. In a similar way we suggest to split the second 
transition metric in two components that can be used for two 
different transitions on the trellis diagram. In this way, at the 
time t — 2, AAfw is approximated with AMW as follows: 
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This metric increment can be easily split on the two transitions 
at time t = 2, which we name as left and right transitions 
supporting transmission at time t. 



To improve decoding, we also propose another approxi- 
mated metric whose computation requires an interaction with 
the Viterbi algorithm: in other words, with this second ap- 
proach the left part of the metrics at time t — 2 can be 
computed by using the outcomes of the Viterbi algorithm after 
having processed left transitions at time t. If the metric com- 
putation runs in parallel with the Viterbi algorithm, when we 
know the survivor path at each state of the trellis between left 
and right transitions, we can compute the metric increments on 
the right transition using the information on the survivor path 
at each departing state a. The second approximated metric can 
be written as 
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where an additional parameter (3 is included for further opti- 
mization. 

We have compared these two metrics and we had confir- 
mation that the second one performs better, since it exploits 
the additional information provided by the Viterbi algorithm. 
In the next section, we will address the way of generalizing 
this method of evaluating decoding metrics for different codes 
with BPSK modulation and different puncturation matrices. 

III. Decoding metric computation 

When the transmitted super-symbol carries coded bits from 
two different transitions due to puncturing, the approximated 
metric increment that exploits the information on the survivor 
path has proved to give the best performance. However, up 
to now, it can only be computed for a puncturation pattern 
similar to that given in Q. The next step is to generalize the 
metric computation for the generic puncturation matrix and for 
a generic number of transmitting antennas. We will do this in 
two steps: first we will find a generalization for the case with 
two transmitting antennas (N — 2) and a generic puncturation 
matrix and then we will remove this hypothesis and consider 
the case with N > 2. 

We first considered a two columns puncturation matrix as 
follows: 

In this case we have three symbols (5 = 3) carrying coded 
bits belonging to different transitions of the trellis diagram. As 
the decoder computes three metrics using the three received 
symbols, we must split these three metrics into four trellis 
transition. The relationship between the received symbols r' 4 ' 
and the coded bits is shown in Fig[3] 

Since carries bits from the second and third transitions 
on the trellis we will use it to compute the metric increments 
for the left and right transitions at time t = 2. Since also 
r^ 3 ) carries bits from the third transition we will use both 
A 2 ^ and A 3 ^ to compute the metric for the third transition 
(which is also the left transition at time t = 3). On the left 

(3) 

metric computed at time t = 3 the symbol c\ is unknown; 



therefore this left metric is evaluated with the value of c 



(3) 
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Fig. 3. Relationship between received symbols and coded bits of P 2 -STC 
obtained from a R=l/2 convolutional code with N=2, M=l 

and puncturation matrix as in (9). 



that minimizes it. On the right metric computed at time t = 2 
the information from the Viterbi algorithm can be exploited 

(2) 

to derive the value of c 2 for the computation of the metric. 
Each metric increment coming from received symbols carrying 
bits belonging to two trellis transitions can be split into two 
components, a left metric AmjJ' and a right metric Aifi^' : 
the former will be approximately computed by considering the 
coded bits of the right transition at time t as unknown, whereas 
the latter will be computed by considering the coded bits on 
the left transition at time t as those of the survivor path at the 
trellis state er (t) . That is, for t = 2, 3, 4: 



AM''' = Ami'' + Am 



(9) 



where 



where is the number of not erased coded bits in the label 
of the left transition at time t and tir is the number of not 
erased coded bits in the label of the rigth transition at time t. 
For example, if we refer to matrices in (20), we suggest using 
= 1/3 in the approximated metric at time t — 2. 

A. A simplified metric 

We propose here another possible approach for the general- 
ization of the metric computation, which is sligthly simplified 
with respect to the one already presented. We name this 
approximated metric as Type-2 metric to distinguish it from 
the former metric named Type-1 metric. With this second 
approach, if we refer to the example in Fig(3]as before, we will 
exploit the information of the survivor path at each state only 
in the last transition with an erased bit, i.e the right-metric 
increment is used only at the fifth transition. As we want our 
decoder to be an approximation of an ML decoder, we have 
to weight each component of each metric increment with an 
appropriate parameter: the first 8 components (left metrics) 
and the last component have to be multiplied by parameters 
whose sum must equal to 6 (equal to 3 in the example). Instead 
of using 8+1 different parameters we choose to use only 
two different weights: the first one, cj a , is used for the left 
metrics, whereas the second one, uib, is used for the last right 
metric. Once again we include a parameter to balance the 
two different weights. 
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In the decoding algorithm for this example, the second 

(2) 

transition of the trellis will be labelled with Arh L , the 

(2) (3) 

third transition with Am K R ' +Am L , the fourth transition with 
Am^'+Am^ and the last transition with Am|Jl 

We have investigated the optimum value of for different 
values of 8 and for a rate 1/2 code with puncturation matrices 
with 2 zeros (as in eq. (9)). We found that the optimum value 
of is a decreasing function of 8. The reason of this behaviour 
is that the information held by the survivor path at each state 
exploited in the right metric will be more and more unreliable 
as the 8 increases. For 6 = 1 this value approaches 0.5. 

The generalization of this approach to the N > 2 case can 
be easily done if we assume that the puncturation pattern does 
not allow the transmission of super-symbols carrying coded 
bits from more than two trellis transition. Also in this case 
each approximated metric increment can be split into a left 
and right components. These two component may be related 
to more than one coded bits and the value of f3 should be 
suitably chosed and optimized for each puncturation pattern, 
for each value of N and even for each time t. However, we can 
give the following rules of thumb to compute the parameter 
(3 in the approximated metric increment at time t: 
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u b = f3- 
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5 + 0(1-8) 

j5 

'6 + 13(1-6) 

With this choice of oj a and u> b , if = then uj a = 1 and 
ujb = 0, i.e. we do not trust the right metric using survivor 
paths information. If = 1 then uj a = and cob = 1, i.e., we 
only exploit the last right metric If — 0.5 then uj a = ujb, i-e. 
we weight the metrics uniformly. For our decoding rule to be 
an approximation of the ML decoder we require that: 



5uj a + u>b = 6 



(15) 



which is clearly satisfied by our construction of parameters uj a 
and u>b- 

If we refer to the same example as before, the expression 
for the approximated metric increments at time t = 2,3,4 
becomes: 



^(4) 
R 



(16) 



t=2,3,4 



(=2,3,4 



where 
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With this approach we have split 6 metric increments 
corresponding to the 6 symbols received at t = 2,3,4 into 



<5+l metric increments. As in the previous case the optimum 
values of j3 is a decreasing function of 6 since the state is more 
and more untrustworthy as 8 increases. A reasonable value of 
(3 is between 0.5 and 1. 

In rate-adaptive applications this allows to change simple 
parameters for metric computation in the Viterbi decoder 
to decode very different puncturation matrices, keeping the 
overall trellis structure of the decoder unchanged. 

IV. Search for good puncturation matrices 

We address now the issue of finding good puncturation 
matrices that allow to design the wanted rate R for the P 2 - 
STC starting from a convolutional mother code of rate 1/N 
for BPSK modulation (this is the case addressed in this paper: 
we can easily extend the approach for any pragmatic STC, i.e. 
as an example by considering a mother code with rate 2 / (2 N) 
or 1/(2N) for QPSK modulation). 

As a first step, we investigate the behavior of simple punc- 
turation patterns of N zeros to be used as building blocks for 
good puncturation matrices with N rows (N is also the number 
of transmitting antennas). We conjecture that the position of 
a basic pattern inside a puncturation matrix do not affect the 
performance (or its effect is negligible). 

A simple observation will guide our search for good pat- 
terns: if the number of super-symbols which carry coded 
bits of different transitions of the trellis diagram increases, 
the performance of the decoder get worse because of the 
increasing for the number of approximated metric increments. 
On the other hand, if we limit this number to zero by 
considering a puncturing pattern of N zeros that erases all the 
bits in a single transition (in this case the entire super-symbol 
is erased), the performance may be degradated because all the 
erased bits are concentrated in a single position. We consider 
basic patterns obtained by placing N zeros on one or more 
subsequent columns (in the latter case d is the number of 
columns). The simulations have proved that putting all zeros in 
the same column is not the best solution. We made simulations 
comparing the performance obtained by each different pattern 
when parameter (3 is fixed to the optimum value for each case. 

For N = 2 we obtained as best puncturing patterns the 
following: 

'=(::: J ?:::)'= (:!J::) <"> 

For N = 3 all the best patterns have one zero on the first 
column and two on the second column, e.g. 
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(20) 



As a second step in our search, we try to use this basic 
patterns to build up more complex puncturation matrices that 
allow to reach higher code rates. In the search for good 
puncturation matrices we put the zeros in a way to satisfy the 
rate-compatibility rule. In table [IV] we show a family of rate- 
compatible puncturation matrices for the N = 2 case with 
a puncturation matrix with p = 10. In table [II] instead, we 
addressed the case with three transmitting antennas (N = 3) 
and a puncturation matrix with p = 10. The final rate R of the 



code will be p/[(m — p)N], where m is the number of basic 
puncturing patterns of N bits in the matrix. 

TABLE I 



rate 


rate-compatible 


R=5/9 


1101111111 
1011111111 


R=5/8 


1101011111 
1010111111 


R=5/7 


1101010111 
1010101111 


R=5/6 


1101010101 
1010101011 


TABLE II 


rate 


rate-compatible 


R=10/27 


110 1111111 
110 1111111 
1011111111 


R=10/24 


1010111111 
1010111111 
10 1111111 


R=10/21 


1010101111 
1010101111 
10 10 11111 


R=10/18 


1010101011 
1010101011 
0101010111 


R=10/15 


1010101010 
1010101010 
0101010101 
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Fig. 4. Effects of different basic puncturation patterns of 2 zeros on the 
performance of the P 2 -STC scheme with N = 2. The mother code is( 133, 17 1) 
R=l/2, with a puncturing period p = 10 



V. Numerical results 

The Figj4] shows the sensitivity of the performance of a 
P 2 -STC with the same rate to different puncturation patterns 
of 2 zeros and to the two types of approximated metrics 
discussed before. We can easily note that Type-I metric is 
better than Type-II metric as expected, and in general, as 
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Fig. 6. Performance over a BFC of a R = 5/8 P 2 -STC obtained from 
rate 1/2 convolutional code (133,171) with N=2,M=1. a)L=10, b)L=8, c)L=6, 
d)L=4, e)L=2 



appropriate choice of the metric increments a single Viterbi 
algorithm on the same trellis diagram can still be used for 
decoding different rate-compatible codes. We have proposed 
two different approximated metrics and explained how to 
construct rate-compatible puncturation matrices that give good 
performances. In the results we have shown that with this 
scheme we can easily obtain codes with rates higher than 1/N 
which can approach or achieve the maximum allowed diversity 
degree on both quasi-static channel and block-fading channel. 
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already commented in previous sections, the performance 
degrades as i5 (distance between the 2 zeros) increases. The 
best case is 5 = 1, even though for Type-II metric performance 
changes sligthly. 

In Fig|5]we show that it is possible to achieve the maximum 
diversity allowed on a quasi static channel with a P 2 -STC 
with 3 transmitting antennas: for 1/3 < R < 2/3, the 
maximum diversity degree is 2 and can be achieved with all 
the codes with rate larger than 10/21. It is evident that there is 
a threashold effect around rate 0.5: if we puncture too much, 
we lose 1 degree of diversity. 

In Fig|6] we show the performance over a BFC of a rate 
5/8 P 2 -STC obtained from the well-known rate 1/2 mother 
code (133,171) with a puncturing pattern like the one in ([8]) 
for varying fading levels per codeword. This code, which can 
not achieve a diversity degree larger than 1 on a quasi-static 
channel, is able to capture a significantly large diversity degree 
according to the limits given by (1). 

VI. Conclusions 

In this paper, we have shown that it is possible to use 
punctured convolutional codes to build pragmatic space-codes 
over a block-fading channel. We have also shown that by an 



